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ABSTRACT 


The apparent diurnal Martian surface pressure variation, as 
deduced from radio occultation experiments, is discussed and 
explained as possibly arising from the effect of a low altitude 
electron layer. Possible source and loss mechanisms for the low 
altitude electron layer are presented and discussed. Time-dependent 
differential equations describing the electron layer are derived and 
then integrated to Investigate the electron distribution resulting 
from the several processes that might occur in the atmosphere. It 
is concluded that the source mechanism is the sublimation of alkali 
atoms from a permanent dust layer (a dust layer of 0.2 micron particles 

V 

_ o 

of density 9 cm J is sufficient), and that the dominant loss process 
must' involve CO 2 clustering to the alkali atoms. Using these processes, 
an electron layer is developed which would explain the apparent diurnal 
surface pressure. 
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CHAPTER I. Introduction 

Since 1965* there have been four U.S. planetary space probes used 
in the exploration of the planet Mars. These spacecraft were used to 
gather data concerning the physical properties of the planet, and one 
of the pr imar y areas of research was the determination of the physical 
properties of the Martian atmosphere. These properties were examined 
by techniques of visual imagery (refs. 1,2,3), infrared spectroscopy 
(refs. ^*,5,6), infrared radiametry (refs. 7,8) ultraviolet spectrome- 
try (refs. 9,10,11,12), and radio occultation (refs. 13, l 1 *, 15, 16,17, 

18,19). 

Analyses of radio occultation data from the U.S. Mariner space- 
craft have produced a considerable amount of information on both the 
neutral atmosphere and ionosphere of Mars. Data from a single occul- 
tation entry and exit were obtained from each of the flyby missions of 
Mariners 4 ( 1965 ), 6 and 7 ( 1969 ), and data from a large number of oc- 
cul tat ions of the Mariner 9 (1971) orbiter have greatly extended the 
latitudinal, longitudinal, and diurnal coverage of the planet. 

On the basis of the currently available analysis of the Mariner 
radio occultation data (refs. 16,17), there are some indications of 
apparently systematically higher surface pressures on the nightside of 
Mars than on the dayside. This difference is approximately 1 mb or 
roughly 20 percent of the total pressure. These indications of an 
apparent diurnal variation in surface pressure might be explained by 
either a variation in topography at the occultation location, or by 
some sort of diurnal process occurring in the atmosphere. This paper 



will discuss one atmospheric diurnal process that might explain the 
radio occultation results. 

In a radio occultation experiment, as the spacecraft passes be- 
hind the planet the atmosphere acts as a lens to the radio ray and 
changes the apparent motion of the spacecraft. The changes in motion 
of the spacecraft are measured by the change in frequency of the re- 
ceived radio ray (Doppler signal). After subtracting the predicted 
change in spacecraft motion from the measured Doppler change, the dif- 
ference or residual is used to deduce the atmospheric density. It has 
been shown that there are interpretation ambiguities inherent in the 
reconstruction of atmospheric properties from single frequency radio 
occultation data (ref. 20). This ambiguity arises from the fact that 
the deduced refract ivity can be composed of both a positive component 
(due to the neutral atmosphere which causes an apparent motion away 
from the observer) and a negative component (due to the presence of 
free electrons which cause an apparent motion towards the observer). 
Thus, a single frequency occultation measurement can be represented by 
one equation with two unknowns which cannot be solved for the effects 
of the two components separately. Because of this interpretation 
ambiguity, numerous atmospheric models (representing different atmos- 
pheric states) can be developed that will produce the refractivity 
profile obtained from any single frequency radio occultation experi- 
ment. One atmospheric model that accounts for the apparent pressure 
disparity includes a time-dependent low altitude electron layer in the 
dayside atmosphere (ref. 26 ). 
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The properties of a planetary atmosphere and ionosphere vary with 
altitude, latitude, longitude, and the local time of day. A low-alti- 
tude electron layer will possess similar variations. This paper will 
discuss an atmospheric model that includes an electron layer, and will 
develop equations to describe the altitude and time variations of the 
layer. To do this, the simplified (no transport) time-dependent elec- 
tron density continuity equation is normalized and non-dimensionalized 
for ease of handling. The resulting first order, non-linear, differ- 
ential equation is integrated to give diurnal electron density profiles 
as a function of various absorption and recombination coefficients. 

The development of these models is somewhat the reverse of that associ- 
ated with the classical Chapman electron layer theory. In that theory, 
the _ atmospheric properties are known and the resulting electron layer 
is determined. In the present case, the magnitude and functional form 
of the electron layer at a specific geographic location and local time 
is specified by the difference between the refractivity resulting from 
the neutral atmosphere and the refractivity as measured by radio oc- 
cultation. Knowing the electron distribution then allows the distri- 
bution of the ionizable constituent to be determined. The recombina- 
tion and the absorbtion coefficients in the atmospheric model are 
empirically adjusted until an acceptable ionizing constituent model 
which produces the specified electron distribution results. 

The next section of this paper will provide a brief description of 
the radio occult at ion experiment and the methods by which the physical 
properties of the Martian atmosphere are inferred from the occultation 



k 


experiments. Following that is a section that discusses the diurnal 
pressure variation as it appears in the existing data, and another 
section that describes the use of that variation to estimate the 
functional form of the proposed electron layer. There are then sev- 
eral sections that discuss the electron layer — its description, 
origin, and methods for modeling it. The remainder of the paper sum- 
marizes the derivation of the equations used in the model, and presents 
the results of the numerical integration of those equations. On the 
basis of these results, conclusions concerning the electron layer and 


the inferred source distribution are drawn. 



5 


CHAPTER II. Radio Occult at ion and the Martian Atmosphere 

The radio occultation occurs as a spacecraft passes behind a 
planet as viewed from the tracking station. The radio signal between 
the spacecraft and tracking station is changed by passing through the 
atmosphere of the planet being studied. By making certain assumptions 
about the atmosphere, the position of the spacecraft, the propagation 
paths, and the change in the received radio frequency (Doppler) as the 
radio ray passes through the atmosphere can be related to the index of 
refraction of the atmosphere. Since the index of refraction ( p ) of 
most gases is numerically very close to unity, the quantity generally 
used in describing the atmosphere is the refractivity unit. The re- 
fractivity is related to the index of refraction as: 

N = (p - 1) * 10 6 . 


Figure 1 is a sketch of a typical profile of refractivity versus 
altitude for the atmosphere of Mars. Plots such as Figure 1 are the 
usual output of radio occultation experiments, and it is these plots 
that form the basis for the inference of atmospheric properties. The 
negative refractivity peak which occurs at 135 km on the dayside pro- 
file is attributed to free electrons resulting from photoionization 
processes involving CO^ in the Martian atmosphere. If there are at 
most a very small number of low mass ions present (a reasonable as- 
sumption for the atmosphere of Mars), the number of free electrons is 
simply related to the negative refractivity by a linear equation. 
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The lower positive peak is assumed, by most authors, to be caused 
only by the neutral atmosphere. The relationship between the refr ac- 
tivity and the properties of the neutral atmosphere is more complex 
than that for electrons since there can be several gases which contrib- 
ute to the refract ivity and only one equation relating them. Therefore, 
the exact gas composition must either be known or assumed, and the re- 
fractivity profile is used to determine the molecular number density. 

The molecular number density is used together with the hydrostatic equa- 
tion and the perfect gas law to arrive at the temperature distribution 
in the atmosphere. The pressure variation in the atmosphere can be ob- 
tained by either integrating the number density in the hydrostatic 
equation or by using the temperature distribution in the perfect gas 
law., 

The peak dayside electron density, as deduced from the Mariner 
radio occultation, occurs at about 135 km with a magnitude of approxi- 
mately 1.6 x 10 5 electrons per cm^. The deduced surface temperatures 
vary from l4l°K to 272°K with about 2Uo°K being the average dayside 
surface temperature and about l60°K being the average nightside surface 
temperature. The deduced lapse rate in the atmosphere varies from 0 
to + 3.8°K/km, with a large number of measurements implying a near 
isothermal, atmosphere, and almost all measurements having a lapse rate 
of less than half the theoretical adiabatic lapse rate (5° K/km) (ref. 
19). (These small temperature gradients are the basis of using, in the 
following sections, the approximation that the atmosphere is isothermal 
very near the surface. The error in pressure using an isothermal 
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approximation over a height of 10 km vill he no larger than about 3 
percent at 10 km. ) The deduced Martian surface pressures vary from 
2.5 mb to 10.8 mb, with the average surface pressure being from ^.7 to 
6.1 mb, depending on what surface is taken to be representative of the 
mean planet surface. 
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CHAPTER III. Description of the Apparent Diurnal Variation in the 
Surface Pressure 

On the basis of the currently available analysis of the Mariner 
radio occultation data (refs. 16,17,18)* there are indications of ap- 
parent systematically higher surface pressures on the nightside of 
Mars than on the dayside. It is difficult to examine the published 
occultation surface pressures as a group, since the data occurs at 
different altitudes, places, times of day, and time of year. To com- 
pare the data as a group, the effects of different altitudes and tem- 
peratures in the atmosphere are taken into account, and any of the 
previously mentioned effects (including season and meteorology) will be 
assumed to cancel out when group averages are taken. 

- A reasonable approximation for a Martian equipotential surface can 
be obtained from the dynamical flattening (ref. 21). An approximation 

to the equipotential radius as a function of latitude is R 2 R [l - 

eq 

p 

f sin (LAT)], where is the equatorial radius, f the flattening co- 
efficient, and LAT is the latitude of the point in question. From 

Mariner 9 orbital analysis (refs. 22,23) the flattening coefficient, 

-3 

f, was measured to be 5.25 x 10 , and this value will be used in the 

following analysis. 

To allow for the different values of deduced surface temperatures 
and altitudes associated with each deduced surface pressure, the de- 
duced surface pressures will be compared by two methods. The first 
method is to scale each surface pressure data to a pressure altitude. 
The pressure altitude is the altitude in a reference atmosphere at 
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which a given pressure level occurs. To illustrate this, assume that 
the Martian equipotential surface has the same dayside and nightside 


pressure (P q ) at the zero altitude or base surface. 


then for a deduced 
RT. 


pressure, P^, temperature , (associated scale height, EL 

altitude , Z . , 

* 1* 



) , and 


P. = P exp 
i o r 



The reference atmosphere is taken to have the seme equipotential sur- 
face, the same zero altitude pressure, P q , a constant temperature, T , 
(and, therefore, scale height, H q ) over the entire reference surface. 
Then, for the same deduced pressure, P^ there is a pressure altitude, 
Z q £., such that 

P i = P o ex P {- sfj - 

and therefore, the pressure altitude associated with each P^, T^, and 
Z^, deduced from a radio occultation experiment, is defined as 

H 

„ _ o „ 


Now, if the base surface pressure is the same for both the day- 
side and nightside surface, then the In of the pressure data plotted 
against the pressure altitude should all lie on a straight line (with- 
in the measurement accuracy and the assumption of an isothermal 
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atmosphere over the heights in question). 

Figure 2 is a plot of In of pressure versus pressure altitude 
above a geopotential ( R e(Juator = 3390 km, f = 5.25 x 10 - ^) for a set 
of night (exit) and day (entrance) deduced surface pressures. A ref- 
erence temperature of 2Uo°K was used and therefore H q = 12 km and the 
deduced surface temperature at each data point was used for IK to 
calculate the pressure altitude. As can he seen from Figure 2, the 
nightside (exit) pressure measurements cluster above the line through 
the dayside pressures, implying the assumption of an equal base sur- 
face pressure, for both the dayside and nightside, was invalid. 

The second way to compare the data is to examine the base surface 
pressure P Q by the relation 

p o = r. exp | + , 

where P^, and IK are the surface pressure, altitude from the geo- 
potential, and surface temperature deduced from the radio occultation 
experiment. Figure 3 is a plot of the calculated base surface pres- 
sure versus relative time of day for a set of day (entry) and night 
(exit) measurements referenced to a geopotential. Again the differ- 
ence between day and night pressure measurements is evident and is 
approximately 1 mb. 

Since it is considered unlikely that a day-night pressure differ- 
ential of such magnitude [20% of the total pressure) could actually be 
sustained for any appreciable time, some process or combination of 
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out that if there are free electrons below 50 km on the dayside, 
these electrons would contribute a negative component to the refrac- 
tivity in this region. The observed refractivity would then be a com- 
posite of the negative refractivity due to the free electrons and the 
positive refractivity due to the neutral atmosphere. The observed 
composite refractivity would then be smaller than that of the actual 
neutral component of refractivity, and therefore, the pressure deduced 
from the observed refractivity would be less than the actual atmos- 
pheric pressure. 
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CHAPTER IV. Estimation of the Form of the Electron Layer 

If it is assumed, that the dayside positive portion of the lower re- 
fractivity profile, as shown in Figure 3, is a composite profile of the 
neutral atmosphere plus an electron layer, then it can he shown that 
real part of the index of refraction is (see Appendix A): 

K Mi) 


where P c is the measured composite real index of refraction, P n is 

the real index of refraction of the neutral atmosphere, and the re- 
maining terms are the contribution of the electron layer (where n is 

the number of electrons, e~ and m are the charge and mass of an 

electron, e Q is the dielectric constant of free space, and f is 

the frequency of the probing radio signal). 

Rearranging Equation (l) and introducing refractivity N = (p - l) 




M2) 


where N g is the refractivity of the electron distribution, N c is 
the measured refractivity, and K q is the refractivity of the neutral 
atmosphere. Equation (2) can be approximated to: (see Appendix A) 




M3) 



lU 


Assuming that the atmosphere is in hydrostatic equilibrium and 
that the neutral refractivity is proportional to the neutral density 
(a valid assumption for the low gas density on Mars), then 

«*/>[- z ^ 4 j u-(U) 

where Z is the altitude above the geopotential and is the scale 

height of the neutral atmosphere 

H. = i / '»/' '- ,5 > 

By assuming that the deduced nightside surface pressure is more indica- 
tive of the actual surface pressure, and by allowing various formula- 
tions for H n , i.e., isothermal, constant temperature gradients, etc., 
many models of electron distributions can be derived. Figure *» is a 
sketch of electron number density versus altitude for several derived 
electron distribution models. The upper solid curve is the ionosphere 
electron layer that is deduced by radio occult at ion for most dayside 
occultations. The lower curves are the calculated electron distribu- 
tion models which would essentially equalize the surface pressure, and 
would still yield the same refractivity profile as the refractivity 
profile observed by the radio occultation experiment. As can be seen 
from Figure U, these calculated layers all have peak densities of 

i 3 

6 to 8 x 1CT electrons per cnr and the peak density occurs at alti- 


tudes below 15 km. It should also be noted that the calculated 
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electron density profiles are only rough estimates, since the data 
from which they are calculated consist primarily of published numbers 
of the pressure at the point of radio occultation and the measured re- 
fract ivity profiles in the atmosphere (of which very few have been 
published). To compound the problem there are large differences in 
the surface pressures, deduced from the same radio occultation data, 
by different authors. "The difference ... is due to the method of 
removing from the Doppler residuals, drifts, and oscillation, whose 
presence is not understood by the Mariner experimenters." (See refer- 
ence 27.) 



CHAPTER V, Bouyce of the proposed fflectron Layer 

In the past, Martian ionospheric models have been extremely con- 
troversial, and a considerable amount of research has teen done on the 
upper ionosphere of Mars (refs. 13,28,29,30,31). Although the exist- 
ence of low electron layers on Mars has been conjectured (ref. 31) and 

to «3 s 

electron density of the order of 10 cm peaked at 65 km have been 
measured (ref. 32), little research has been reported on the lower 
ionosphere, Whitten, et al, (refs, 33,3to) have studied the lower 
Martian ionosphere by studying the influx of solar protons and cosmic 
rays on a C0 2 atmosphere. They concluded that for a quiet sun the 
lower ionosphere below 65 km is formed predominately by galactic cos- 
mic rays, and is a permanent layer with little diurnal variation 
peaked at approximately 25 to 35 km with a maximum density of the 

o o 

order of 10 J electrons cm • 

Since all of the efficient channels of electron production in- 
volving photochemical processes in gaseous carbon dioxide in the 
Martian atmosphere have been utilised in theoretical modeling of the 
ionosphere, the postulated electron layer must come from some addi- 
tional process in the atmosphere. Also, since the effect of the 
electron layer is diurnal, the electron density should exhibit a 
diurnal variation. Assuming that this diurnal variation is driven by 
sunlight the source must be able to liberate electrons when exposed to 
solar ultraviolet radiation. It is known that COg gaa absorbs sig- 
nificantly in the 100 ~ i?00£ 6 region (ref. 25) and CD, ; -r. -*> 

w- - 

plstely transparent from about 173>0&° at least 21CCA 0 I [ref. , 


IT 


The Martian atmosphere has very little ozone and therefore there is an 
ultraviolet window in the 1750A 0 to 2100A° region (ref. 35). Hence, 
light in this wavelength interval can reach the Martian surface. It 
follows then that the electron source material must release electrons 
when exposed to light in the l800A° to 2100A° region ( 6.9 ev to 5*9 ev). 
Based on ionization potentials and abundance criteria, the list of 
source candidates has been reduced to: potassium (ionization poten- 

tial = U.3 eV), sodium (5-12 eV), barium (5.19 eV), lithium (5.36 eV), 
aluminum (5.96 eV), and calcium (6.09 eV). All of these materials are 
cosmically abundant and found in both planetary crustal material and 
meteoric dust. 

This electron source material, alkali metals in the atmosphere, 
must come from either the top of the atmosphere (meteoric ablation), 
the bottom of the atmosphere (crustal material raised into the atmo- 
sphere by the wind), or a combination of the two methods. It is pro- 
posed that the alkali metals sublimate from aerosols or "dust" par- 
ticles that have arrived in the atmosphere by one of the methods just 
mentioned. 
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CHAPTER VI, PuBt and Aerosols in the Martian Atmosphere 

Dust sad aerosol particles in the Martian atmosphere have been ob- 
served and measured by many authors (refs. 37 thru 1*3) . Thin detached 
haze layers were distinguished in Mariner photographs at altitudes of 
5 to 1*5 km (ref. 1*1* ). The presence of dust in the atmosphere was 
found to be able to account for the fact that measured atmospheric 
temperatures were considerably warmer than the theoretical tempera- 
tures calculated, assisting radiative equilibrium in the atmosphere 
(ref. 1 * 5 ). The presence of a low altitude dust layer would also ac- 
count for the observed photolytic stability of the Martian CO^ 
atmosphere (ref. 1*6) and in addition this dust brings into agreement 
the observed and computed profiles of the minor constituents 0, C> 2 , 

00 and 0^ (ref. 1*6). In fact some authors have concluded that there 
are aerosols, that absorb solar radiation, permanently present in the 
Martian atmosphere (ref. 1*7). This dust would most likely be composed 
of oxides of Si, Ti, Fe, Mn, Mg, Ca, and Ha (ref. 1*8). The dust in the 
Martian atmosphere has a measured SiO^ content of about 60 percent 
(ref. 39); and since this is a slightly enhanced silicon content, the 
alkaline content will also be enhanced (ref. 1*9). Therefore, if a 

.2ym particle lost 1/10,000 th of its weight in sublimated alkali 

k 

atoms, it would release on the order of 2 * 10 alkali atoms into the 
atmosphere. 

The dust could arrive in the atmosphere in either or both of two 
ways. First, the dust may be the result of the surface material having 
been broken up into fine powder by the extreme heating and cooling 
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differences experienced over a Martian day. This powdered surface ma- 
terial could easily he blown into the atmosphere by thermal winds (ref. 
50,51) or localized dust storms which occur rather frequently (ref. 

52). Since .2y particles would remain in the atmosphere on the order 
of two years, it is concluded that the atmosphere may contain such 
particles most of the time (ref. 53). Second, the dust could result 
from the ablation of meteoroids which enter the atmosphere. This idea 
has been advanced as the source of sodium layers in the earth (ref. 5*0 
and since Mars is closer to the asteroid belt than earth, the effect 
may be even more pronounced for Mars. 
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CHAPTER VII. Modeling of the Alkali Layer 

Aerosol layers have been detected in the earth's lower atmosphere 
(refs. 55*56,57). Alkali metal layers have also been measured in the 
earth's atmosphere (refs. 5^,58,59*60,61 ,62) . The sodium layer in the 
earth's atmosphere has been postulated to come from an aerosol or "dust 
bank layer" (refs. 62,63,61+), or from meteoric material (refs. 5^,65* 
66). The actual photochemical process and distributions in the earth’s 
atmosphere are not too well understood or modeled, primarily because of 
the lack of quantitative atmospheric data, lack of laboratory measure- 
ments of the physical properties of metal ions and their oxides, and 
lack of knowledge as to the effects of competing reactions such as 
C0 2 clustering (ref. 67 ) and aerosol ion pair annihilation (ref. 68). 

- It will be assumed that the final Martian alkali number distribu- 
tion follows roughly the same shape as that found for earth's sodium 
distribution (such as shown in reference 65 ) and that this distribu- 
tion can be adequately approximated by the equation: 

n a (2,0) = K e *p /- - 0 e. 7-d) 

where n^(Z,0) is the alkali number density as a function of altitude 
(Z), at time, t = 0. A can be thought of as a "scale height" for 
the rate of fall of the upper side of the distribution and B is a 
dimensionless constant which controls the shape of the lower side of 
the distribution. As a starting basis, A was given the value of 6 km 
which is approximately one-half the neutral scale height as observed in 
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the earth's sodium layer (ref. 6l). The altitude of peak density of the 
alkali distribution is given by 

£o*„ t-(2) 

Thus, the value of B was chosen to be 6 in order that the distribu- 
tion peak occur at an altitude below 15 km (for the numbers chosen the 
peak occurs at. 10.75 km). 



22 


CHAPTER VIII. Derivation of the Time-Dependent Electron Distribution 
The continuity equation relates the change in electron density per 
unit time to the production, loss, and the divergence, and is given by: 

tv = p-a - v-c«.y; 


where is the electron density, p is the production rate, L is 

the loss rate, and V is the transport velocity. It will be assumed 
that the transport terms can be neglected, so that the only terms to be 
derived are the production and loss rates. Figure 5 depicts the 
geometry and defines the variables used in the derivation. (As can be 
seen from Figure 5 the atmosphere is assumed to be spherically sym- 
metric.) It is assumed that the absorbtion of ionizing solar radia- 
tion in an element of atmospheric pathlength, dS, is proportional to 

the radiation flux, I, the atmospheric alkali concentration, N , and 

a 

the absorption cross section, Q, and thus can be written as 

JI- 8-(i) 


where the I and Q are wavelength dependent. This wavelength de- 
pendence will be accounted for later in this development. 

Now from Figure 5 


and 


-ft - h)ai* i X * (Rf 2 )*/* X 

~ -pCdi X 


8 -( 2 ) 
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2 

where x is the solar zenith angle, and so that dS - pdsc XdX 
which will be used in Equation 8-{l). 

In addition to loss of radiation due to ionization, there is a 
competing reaction from the absorption by CO^ gas whose absorption 
coefficient is 0. Thus, Equation 8-(l) becomes 

<^1* -I + 8— (3) 

It is assumed that the diurnal variation in the concentration of 
the alkali material is very small (in keeping with the long lifetime 
of the particles); and since their number density is also small, the 
removal of a source atom must be accounted for in some manner. 

. Therefore, let 

8-(io 

Initial calculations showed that n (Z,t) follows the n distribu- 

6 A 

tion almost exactly so that it will be approximated as: 


£<*>**> /~* -6e +J - Qct) fch) 


M5) 
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Combining equations 8-(5) and 8~{3) 


<i£ - -Z [q (n^ch t o) -&*) 


8 -( 6 ) 


now from 7-(l) 


Jl = --L /q Fctt JX 8-(7) 


so that 


J “* Q ft*) J f -jy cs£* 


8 -( 8 ) 


or . 






f^LQ&*)(k-(x(jk))+eti 8 -( 9 ) 


where h « pcscA - R 

Now the rate of production of electrons is 


T) 

Js 


8-(l0) 


where A is the number of electrons released per photon absorbed for 
ionization. Equation 8-(9) then becomes 



25 


*P 3 <>) - &ut) 

Mu) 

or 

p- A T a [q (% 0,6) ^M(k-£kti)^ fA^jocdpj 8-(12) 

where h * pcscA - R 

Since the ionization cross section and the solar flux are both 
wavelength dependent, the probability of ionization can be written as 

f* 

-fr* JQlJX 

where Q is the absorption cross section and I q is the solar flux 

-22 

per unit wavelength. For alkali metals Q varies from about 10 
— 18 —2 

to 10~ cm (refs. 69.70) . Integrating over wavelength, P be- 
comes 

TV# 3x# jcc' 7 i 3 x/A^st* 1 

vhere (A^ X g ) - (1900 A°, 2700 A°), 

and these values, which are roughly the same as those for earth (ref. 
71), are used for QI q in Equation 8- (12). 
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The time dependent term in the production term is the variable x 
(solar zenith tingle). The variable x is related to <p (local time 
measured from noon), by 


CoiJt - S/n D W)+ Cost) Cas(unf)Cos (£ 8-(l3) 

where D is the declination of the sun, LAII is the latitude at 
point P at height Z and 

l - /y// 4 > V , since there are 8 - {lU) 

TT f 

88775 seconds in the Martian solar day. 

The conditions for the sun not to be visible at point P are 

< 1 8-05) 

and if these are met 


Cos k <0 /tup j(/i. 


o k 


V - o 8-0.6) 

There were two types of loss reactions considered. The first is 
radiative recombination in which the ionized material recombines di- 
rectly with an electron. This type of reaction is schematically 
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A > A*+C. 


Cp i/ftf •'<>* ) 


8-(l7) 


and 


^ $ "J* Ad t&c&nkisaJ '*** *3 8-(l8) 

If the atmosphere is electrically neutral, this reaction has a loss 
term of 



8 - (19) 


The second type of loss reaction is a more complicated process in 
the atmosphere, such as 

A + +C (pr^doc/t id*) 8 - ( 17 ) 

followed by the two-step process 


A C0J+ + +&>„ (rat. Coefficient *,„) 


8-(20) 


and 


^ ^ (rt-U Coe ff, Clint ) 


8 - ( 21 ) 
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A second reaction of the same type is the process: 

A+t* M -> fl + AA (ret* 

A +e* — > (rat* K?3 ) 


8 -( 22 ) 


8-(23> 


Reaction pair (18-19) are called CO^ clustering (ref, 67 ) and re- 
action pair (20-21) would be using as dust particle as a recombination 
center, such as suggested previously (ref, 68). In both reactions 
(20-21) and (22-23), if the rate constants are such that 


» **3 ^ 


then the loss process can again be characterized by Equation (17), 

4 

The value of a varies greatly, depending on whether the loss 
process is only radiative recombination or some combination of the 


processes discussed above. The range of values used here is from 
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2 * 10~ 12 cm^sec -1 to 10~^ cm^sec 1 (refs. 71 *72,73*710 . 

Finally, the production and loss terms are combined into the 
simplified continuity equation to yield, 


~JT 


~P-L 


8~(2h) 


or combining with Equation 8-(lU) 

where P is defined by Equation 8-(l2), L by Equation 8-(l9), and 
the relation of X to $ by Equation 8-(l3). 

. A computer program was developed to solve Equation 8 - (25). In- 
puts to the program are, height in the atmosphere, latitude, sun de- 
clination angle and constants for the production rate, loss rate, and 
alkali number density. The program starts with zero initial electron 
density and integrates Equation 8-(25) forward in time, using the 
final electron density of the previous day as the starting condition 
for the successive day. 

The differential Equation 8— ( 25 ) is integrated by a fifth-order 
integration subroutine. The classical fourth-order Runge-Kutta 
formula is applied in conjunction with Richardson's Extrapolation to 
the Limit Theory. The subroutine is a variable interval size routine 
in which the interval is varied to meet a specified local relative 
truncation error. A second subroutine is used to compute the integral 
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in the production function (Equation 8-(l2)). The current values of 
al 1 variables, as updated from the differential equation subroutine, 
are used in a ten point Gauss Quadrature integration subroutine. The 
accuracy of these subroutines is far better than is required for the 
problem at this time, since the accuracy of the parameters charac- 
terizing the electron layer are order of magnitude. 

At the end of three Martian days, the results were printed and 
compared. If convergence^ 

AlCtt+oAiAa.*) - /l € Ct* ) _ a 

-- £ / y ^ f ^ 

had occurred, the results were plotted; if not, then the program was 
continued from the last computing point or new intial conditions were 
imposed and the program rerun. 
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CHAPTER IX. Approximate Solution of Some Equations of Interest 

In order to both obtain an estimate of what effect the variation 

of coefficients (such as the absorption and recombination coefficients) 

would do to the electron layer, and to reduce the range and time of 

the computer runs, several approximations to the previously described 

equations were developed. The first was developed to approximate the 

maximum electron density. In order to do this, two approximations 

were made in the electron production term (Equation 8-(l2)); that the 

effect of C0 2 absorption was negligible (a - 10 cm ); and that, 

since the alkali layer has a small number density (characterized by 

the parameter K) and such a small absorption cross section that the 

product QK is sufficiently small enough to drop the exponential time 

dependence term in Equation 8-(l2). Since the maximum electron densi- 

dN 

ty will occur when —rr =0 (if reached), then 




so that if 



hi!! ~ § [ i i + - ij 


9— (1 ) 


9-(2) 


9-(3) 
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From Equation 9~(3) (under the same assumptions as Equation 9-(l)) 
can also he found the required alkali concentration for a given m ax imum 
electron density, i.e.. 


n A (» ~ 




Rn, 


P U) 


R 


9-W 


The second equation is the decline in the electron number density 
after the sun sets. The differential equation is 


d At 

TT 



9-(5) 


the solution of which is 

At*, 

; , 9-(6) 

or rewritten 

- /1«C b,6) 

/?«. -£■>! 

Equations 9-(U) and 9-(T) will he used later to narrow the range 
of acceptable coefficients. The final results to be presented, how- 
ever, are based on the numerical solutions to Equation 8-(25). 


/ 


9— (T ) 
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CHAPTER X* Results and Discussion of the Stuc 


All of the computer simulations were run at a latitude of &nd 
a sun declination angle of 0° (corresponding to summer or winter 
solstice). On the basis of the considerations discussed in Section 

VIII, the range of the absorption cross section was 10 ^ to 10 ^ 

2 -18 
cm and the range of the recombination coefficients was from 2 x 10 

“6 3 *”1 

to 10 cm'sec*’ . 


Table I and Table II list the coefficients for each model and the 


electron densities for each model at eight altitudes in the atmosphere. 
In Table I there are two electron densities listed for each model at 
each altitude. The upper entry is the expected peak equilibrium 
electron density as obtained from approximation Equation 9— ( 3 ) . The 
lower entry is the actual maximum electron density as computed from 
the time-dependent equations. Again, in Table II there are two elec- 
tron densities listed for each model at each altitude. In Table II 


the upper number corresponds to the expected minimum electron density 
obtained by using Equation 9-(6), and the lower number is the minimum 
computed electron density as computed from the time-dependent equa- 
tions. Some general observations can be made from these tables. 

These observations will be shown in further detail in the following 
figures. The first is that Equation 9- (3) gives a reasonable approxi- 
mation for the peak value of the electron number density, and, as 
is shown in equation 9-(3)» the larger the value of R, (R = a 1 /QJ Q dX) 9 
the larger the peak electron number density. The second observation is 
that the smaller the value of an e , the less the nighttime decay of 
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the electron density as was predicted by equation 9-(6). 

Figure 6a and 6b are plots of the electron density for Model 0. 

Model 0 has constants of a = 2 x lo" 12 cm 3 sec -1 , Q = 10 -22 cm 2 , and 

5 -3 

K = 9-78 x 10 cm . In Figure 6a are plotted the electron density 
versus hour angle for eight altitudes in the atmosphere. The hour 
angle is related to the time of day by Equation 8-(l3). Examples of 
local times are <p = -90° is approximately sunrise, ij> = 0 is local 
noon, <p = 90° is approximately local sunset, and <p = - l80° is 
midnight. Figure 6b is a plot of electron density versus altitude for 
the four times just illustrated. Model 0 has a large value of R and 
also has a large peak electron density, 2 x 10^ cm -3 . This model has 
the smallest value of and thus the electron density exhibits an 
almost imperceptible nighttime decay. The lack of decay is exhibited 
in Figure 6b, where the profiles for the four times of day all lie on 
each other. The initial condition for Model 0 was zero electrons per 
cm 3 , but it was converging so slowly that the computation was re- 
started using electron densities slightly lower than the values ob- 
tained from Equation 9“(3). This procedure, in effect, moves time for- 
ward in a large step and results in a convergence of An/n between 

- 1 * 

morning terminators of about 10 

Model 1 (a = 2 x 10 12 cm 3 sec -1 , Q = 10 -1 ® cm 2 , and K = 9.78 x 
10^) results are plotted in Figures 7a and 7b. The value of R for 
Model 1 is very large and this should result in a large electron 
density, which did occur and is shown in Figure 7a. 
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Again the value of an g is small and the variation of electron 
density is very small, which results in all four time curves appar- 
ently coinciding in Figure 7b. Model 1 converged in about a week of 
Martian time. 

Model 2 (a = 10" 6 cm 3 sec" 1 , Q = 10" 22 cm 2 , K = 9.78 * 10 5 cm" 3 ) 
is the first model that shows a day /night electron density variation 
of the type needed to explain the diurnal pressure variation. It 
shows a diurnal equilibrium, i.e., the cycle repeats daily, but it does 
not exhibit an equilibrium in the sense that the time derivative of the 
electron density equals zero. This can be seen in Figure 8a, where 
there is a large discontinuity in d^/dt at the terminators, no 
equilibrium concentration is reached, and yet diurnal equilibrium is 
established. Model 2 has a small value of R, which is exhibited by 
the extremely low electron densities achieved (on the order of 1/500 
of the electron densities of Model 0, for example). The value of ane 
is in the midrange of those investigated, resulting in the slow decay 
in the density distribution during the night. This can be seen in 
Figure 8b which is the first model to exhibit diurnal variation of the 
electron profile. Model 2 was slow to converge and convergence was 
difficult to determine because of the small numerical value of elec- 
tron density. 

Model 3 (ot = 10" 6 cm 3 sec" 1 , Q = 10“ l8 cm 2 , K = 9.78 * 10 5 cm" 3 ) 
completes the extremes of Q and a investigated, and also exhibits 
the third type of diurnal variation encountered in the investigation. 
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This diurnal variation is evident in Figure 9a, where there is a very 
sharp rise in electron density at sunrise, which is a result of the 
large value of Q, and equilibrium of dn e /dt being zero is reached 
during the sunlight hours, and then a very sharp decay in the electron 
density at night, which is a result of the large value of an e . From 
Figure 9b, the profile plot for Model 3, it can be seen that for the 
midrange value of R, the peak electron density is not too very large, 
but the large diurnal variation in electron density, which is required 
to explain the pressure differences, is present. Model 3 
reached equilibrium in about 6 days of Martian time. 

To complete this preliminary analysis, the mid-range values of 
—22 2 —8 3 —1 

Q = 10 cm and a = 10 cni sec were investigated. Model 6 

(ot p 10 cm^ sec Q - 10 22 cm 2 , K = 9»?8 x 10^ cm results are 

plotted in Figures 10a and 10b. The value of R is small, which re- 
sults in low electron number density as can be seen in Figure 10a. 

The value an e was also small, which implies slow decay, which is 
exhibited in both Figures 10a and 10b. Model 6 reached equilibrium 
after about 30 diurnal cycles. 

The results of Model 7 (a = 10~^ cm^ sec^ 1 , Q = IQ” 1 ® cm 2 , K = 

5 3 

9.78 x 10 cm ) are plotted in Figures 11a and lib. The recombination 
coefficient a is such that cm is large and thus the distribution 
should exhibit a pronounced nighttime decay, and this decay can be 
seen from Figure 11a. The value of R for Model 7 is also large and 
therefore the peak electron number density is large as shown in Figure 


lib. 
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Examination of the foregoing results indicates that the sinu- 
soidal diurnal variation of the electron density which is observed in 
the earth* s ionosphere is not duplicated by this low-lying electron 
layer on Mars. This is predominately due to the affect of the product 
of the absorption cross section (low) and the source number density 
(very low) being a very small number. When this product is small the 
exponential term in the electron production Equation 9-(3)» which ac- 
counts for the time variation, is very weak. Therefore, even for the 
long slant paths at the terminators, the ionizing radiation on Mars 
has only been slightly attenuated. In fact, the production function, 
for the low-lying electron, can be thought of almost as a light switch 
having two positions — on and off. 

_ The second result of this study came after the conclusion of all 
the actual time dependent calculations. The result was that Equation 
9-(3) , which neglected time dependence, gave a good approximation to 
the peak electron number density (6 % average error). Thus as more and 
better data about the Martian atmosphere become available, equation 
9— ( 3 ) can be used to give a quick approximation as to the effect of 
the low-lying electron layer on the atmospheric properties deduced by 
radio occult at ion and vice versa. 

The third result was that the larger the value of an the 

© 

faster the nighttime decay of the electron density. This would be ex- 
pected from an examination of Equation 9-(6). The general trend was 

that for values of an < 10 ^ sec 1 the distribution showed no 

e 

diurnal variation; for values 10~^ sec* -1 > an > 10~^ sec -1 , the 
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distribution exhibited a sawtooth variation; and for an g > 10 
sec \ a sharp rise and fall or square wave. 

Combining the results of the foregoing study with the measured 
pressure data should allow the possible range of the various coeffi- 
cient to be narrowed. From the discussion of pressure data is 

Section IV, it is clear that the peak electron density must be on the 
1+ 3 

order of 6 * 10 electrons per cm , and that this number must de- 
crease at night to the point of being undetectable. To determine the 
recombination coefficient Equation 9-(7) is used 


/lcC*j*6) i 
0 < - ~ 

Ac <*,&)£* -kU 


9— (7 ) 


with n g (10.75 Km. sunset) " 6 » 10^ cm {t - t sunset) ~ 2 hours 

and n £ (10.75 Km, t)16» 10^ cm”l which is an electron density low 

enough to be masked by an experiment error of h% at the pressure 

levels in question. For these values, the resulting recombination 

^3 1 

coefficient must be on the order of 2 x 10 - cnr sec~ ; for ease of 
calculation let a *= 10 cm* 3 sec . Then, to determine the required 


alkali source number density. Equation 9-(*0 is rewritten with 
n A (z) = Kf(z) 
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where f(z) = exp { - (z/6 + 6e -Z/ ^ 6 ) } and R = a -1 / QI Q dX. There- 
fore assuming Q = 10 -1 ^ cm 2 , K is of the order of 2.92 * 10^ cm -3 , 
which means a peak density of alkali atoms of 1.79 * 10^ cm -3 . Having 
a 0.2 micron particle sublimate 1/10,000 of its mass in alkali atoms 
would require a dust distribution of only 9 particles per cm 3 in order 
to produce the required alkali atom density. One outcome of this 

study, then, is that the coefficients of a = 10~^ cm 3 sec -1 , Q = 10 -1 ^ 
2 6 -3 

cm and K = 2.92 x 10 cm , which are within the range of laboratory 
measurements of the processes involved, should produce an electron 
density model which will explain the apparent diurnal pressure varia- 
tion. These coefficients are used in Model 8, and the computed elec- 
tron densities are plotted in Figures 12a and 12b. The effect of the 
larger cross section can be seen in Figure 12a as the extremely rapid 
rise in the electron density immediately after sunrise. The electron 

]p n 

density reaches an equilibrium of about 6 x 10 electrons per cm , as 
required, for the entire daylight period. The medium large product of 
an e produces a rapid nighttime decay as is shown in Figures 12a and 
12b. From Figure 12b can be seen the electron density height profiles 
for four different times in the Martian day. The morning profile has 
already nearly 90 percent of the equilibrium density shortly after 
sunrise. The noon and evening profiles coincide at the equilibrium 
densities. Finally the extremely rapid decay can be seen as the mid- 
night profile has fallen to 1/100 of the equilibrium density. 
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XI* Conclusions 

There are some general conclusions about the low-lying electron 
layer distribution that can be inferred from this study. First, be- 
cause the product of the absorption cross section and the number 
density of the ionized constituent is small, the production term in 
the continuity equation has very little daylight time dependence. The 
result of this lack of time dependence is that the sinusoidal time 
variation in electron density that is observed on earth is not present 
in the low altitude electron distribution. Second, the value of 
an e determines the shape of the electron density versus time 
curves; and thus, for electron layers that are of possible importance 
in the explanation of the radio occultation pressure discrepancies, 

the- value of an must be at least equal to or greater than about 
© 

10" U sec" 1 . Third, for an £ > 10 _i+ sec" 1 . Equation 9-(3) gives a 

good approximation to the peak electron density daytime distribution. 

From the preliminary radio occultation data, it has been found 

that the following values for the coefficients used in this study 

would adequately account for a 1 mb difference between day and 

—8 3 —1 — 18 2 

night pressures: a - 10 cnr sec , Q = 10 cm and K = 2.92 x 

IQ? cm“ 3 . 

This value of the absorption coefficient is in the range of that 
measured for the alkali metals. The value of the recombination co- 
efficient 10” cnr sec ) implies that radiative recombination is 
not the dominant process for the loss of electrons. The possible loss 


mechanism that the dust acts as a recombination site, recombination 
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rate unknown, is also not very probable since the dust density is 
numerically so small. Thus, the most probable mechanism for the loss 
of electrons is that of CO^ clustering to the alkali ion (ref. 67), 
which has a very high rate coefficient. 

Finally, because of an analysis of the following factors: a low 

altitude residual dust layer in the Martian atmosphere has been ob- 
served and measured and is theoretically required to explain tempera- 
ture lapse rates; this dust will be composed of alkali and alkali com- 
pounds having low ionization potentials; and the Martian atmosphere has 
low concentrations of 0, 0^, and 0^ allowing solar ultraviolet radia- 
tion to penetrate to the planet's surface, it is concluded that there 
will be a low altitude electron layer. This study has shown that it 

is well within the range of possible conditions in the atmosphere to 

U 

obtain peak electron densities on the order of 6 * 10 electrons cm” , 
which would be required to explain the observed diurnal pressure dif- 
ference. 

XII. Suggestions for Future Research 

Several areas of research are available at the present time uti- 
lizing currently available data. The first is a comprehensive re- 
analysis of present Martian atmospheric pressure data to determine if 
the effects of the low altitude electron layer are of sufficient 
magnitude to be able to more closely define the physical properties of 
the layer. Secondly, studies can be initiated towards a simple 
analysis of meteoric ablation. These studies would give an idea as to 
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how closely an equilibrium distribution of ablated particles matches 
the required aerosol concentration layer. Differences between the 
two, if any, would give a starting point as to the amount of aerosols 
that must be lifted into the Martian atmosphere by meteorological 
phenomena. 

Finally, it is proposed that this low altitude electron layer 

should be searched for and examined by some experiment. The discovery 

and understanding of such a layer would not only be useful in the 

study of the Martian atmosphere but would also aid in the modeling of 

the poorly understood earth’s alkali metal ion region and sporadic E 

layer, since the Martian analysis would be much simpler (because of 

the much simpler Mars atmosphere). For example, several models of the 

-> *> 

earth’s alkali layer involve wind induced V x b shear layers of 
charged dust particles as a source of alkali atoms or involve wind 
induced V x B layers as a removal mechanism for the alkali ion 

—3 

(ref. 75) . Thus, on Mars where the magnetic field is less than 10 
that of the earth, the analysis of the effects of the lack of large 
magnitude V x B forces in the alkali layer formation would have di- 
rect bearing on earth models requiring such forces. 

Two experiments which could be utilized to detect the low alti- 
tude electron layer are a two-frequency radio occultation experiment 
or a twilight glow phenomena experiment. The two- frequency radio 
occultation experiment (currently planned for the Viking orbiter) will 
yield two values for the index of refraction of the Martian atmosphere 
because of the different probing frequencies (see Equation U— (1 ) ) . 



The subtraction of the two indices of refraction and the knowledge of 
the frequencies will yield a close approximation to the electron 
number density (since the index of refraction of most low density 
neutral gases is not very wavelength dependent). The twilight glow 
experiment would be to search for emission lines in the atmospheric 
twilight. The identification of lines allows the delineation of the 
alkali species and the intensity of the line can be related to the 
species number density. Such an experiment could be implemented on an 
orbiting spacecraft and could quite likely be combined with some simi- 
lar type of emission experiment. 
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Appendix A 

This appendix is a short derivation of the steps needed to arrive 

at equations U— 1 and U-3 in section IV. To derive equation W, it 

is assumed at first that in an electrically neutral, non-conducting 

(no free electrons) atmosphere of dielectric constant e , one of 

n 

Maxwells equations can be written as 
A(I) ^7 ~ 

If now the electrically neutral medium contain free electrons 
(density n^) in addition to the neutral part of the medium, then the 
conducting part of the medium can be characterized by a conductivity 
a and Maxwells equation is supplemented by Ohm's law 

— h * * 

a( 2) X r ^ ^ (eT) v 


so that 


where V is the velocity of electron 


A(3) 

or 

A{k) 


= £*ir+ re 


If it is assumed that the right hand side can be expressed as 


9E 

3t 


some composite dielectric coefficient times 


to reduce the 



equation to the simple form of eq. (l) then A(5) V * H = e rr* + <5* B 

n dt 

9E 

= e irr , To determine the conductivity, assume that there are 

C dt 

negligible collisions so that the equation of motion of an electron 
is 


A(6) 

™ jt* 

A(T) 

—y 

now if E = E e 
o 

A(8) 

V = - 


-iurt 


- C 
co m 


from A(2) 


A(9) 




£ * /? t e" S’ 


co m 


so that 


A(10) 


<r - 


6 

co m 


substituting A(10) into A(5) 



A(12) 


to 6 n + 




<L 


u> m 



hG 


or rearranged 


A(13) 



/fc (tT) * 

U >*/> 1 




If the magnetic field is small such that the magnetic permeability 

of the medium is approximately that of free space, then p -4¥ 

o 

is the real index of refraction and 


which is equation U-l. 

To arrive at equation 4-3 in the text, define 


/?«. (grf 


A(15) 



f* 


then, equation A-l4 becomes 


A(l6) 


4 = A.* - X 




Now the index of refraction of just the electrons is 


A(1T) 


/it =(l + */ e /- X 


that 1 + N x 10” 6 = /I - X A/ 1 - i 
e 


so tnat i + w x±u -/ 1 -x CH. 1 -—X 

v o 


A(18) 



A(19) /V* 


substituting A(l6) into A(l9) 


A(20) A/ t X/<> ) z 


A( 21 ) Now u + U - 2 
c n 


A(22) and 


yU. = / -h AJ t/0 


-t 


so that A- 20 becomes 


A(23 ) AJe. X'0~ C £ ~ (/*- AJu jr/<»' 6 - 
/V e 2? A/ c - A/ n 


j (jU-c. "Mn ')(/ie.+ y dlK ) 


/-/!/„* V»' < )&) 


which is equation 4-3 in the text. 
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Table I. Predicted maximum electron densities (upper entry) and computed maximum electron densities 
(lower entry) , for eight altitudes in the Model atmospheres. 
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Table 2. Predicted morning electron densities (upper entry) and computed morning electron densities 
(lower entry) , for eight altitudes in the Model atmospheres. 
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Figure I. Sketch of typical atmospheric refractivity profiles for Mars. 
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Figure 2. Dayside and nightside surface pressures plotted against pressure altitude. 
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Figure 3. Martian base pressures from IV 
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Figure 4. Calculated lower atmosphere electron density profiles. 
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Figure 5. Geometry used in the derivation of the electron 
production equation. 













Figure !0a. Electron 


of Model 6 as a function of hour angle. 


Altitude, km 



Figure 10b. Electron profile for Model 6. 
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Figure lib. Electron profile for Model 7. 




